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Abstract This paper describes an alternative approach for generating point¬ 
ing models for telescopes equipped with serial kinematics, esp. equatorial or 
alt-az mounts. Our model construction does not exploit any assumption for the 
underlying physical constraints of the mount, however, one can assign various 
effects to the respective components of the equations. In order to recover the 
pointing model parameters, classical linear least squares fitting procedures can 
be applied. This parameterization also lacks any kind of parametric singularity. 
We demonstrate the efficiency of this type of model on real measurements with 
meter-class telescopes where the results provide a root mean square accuracy 
of 1.5 — 2 arcseconds. 

Keywords Instrumentation: miscellaneous • Pointing Models • Methods: 
analytical 


1 Introduction 


Pointing models are widely exploited in telescope control in order to correct 
for mechanical and manufacturing imperfections of the various parts of the 
system. The basic concept of a pointing model is to figure out and quantify 
the difference between the targeted and the apparent position of the telescope 
throughout the observations. This quantification can be done using lookup ta¬ 
bles with some grid spacing as well as using analytical approaches where these 
functions are related to the physical behaviour of the telescope mechanics. 

In most of the pointi ng models, either in the case of equatorial 


mounts ( Soillar et al. . 1993), altitude-azimuth mounts ( Granzer et al. . 2012 : 
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Zhang fc W u. 2001) or complex systems of more telescopes (|Gothe et al .1 . 1 201. 'll) 


the raw, uncorrected positions of the various axes (hour axis, declination axis, 
azimuth axis or elevation axis) are read from high precision rotary encoders. 
Such encoders provide a precision within fractions of arcseconds or even hun¬ 
dredths of arcseconds. T here are o ther types of pointing models that involve 


MEMS accelerometers (IMeszaros et al.l . [2014) instead of such encoders. The 


attainable precision is much lower in the case of these integrated accelerome¬ 
ters (currently it is in the range of an arcminute). Despite the lower precision, 
these solutions also have advantages, including the redundant operation (via 
the various constraints between the accelerometer channels), the resistance to 
unintentional tampering and the possibility of installation without modifying 
or altering the drivers. 


The goal of our work presented here is to provide a strictly mathematical 
approach for finding analytical pointing models which provide an RMS ac¬ 
curacy within a few arcseconds. First, in Sec. [2] we detail the mathematical 
formalism used in the derivation of an isotropic pointing model. This model 
is then extended in Sec. [3] with terms that quantify a generalized form of de¬ 
pendence by performing expansions of spherical harmonics. These expansions 
allow the model to take into account various forms of anisotropic behaviour, 
for instance the torques inducted by gravity or misalignments of the gearing. 
In Sec. [5] we report results of test measurements validating these models, as 
well as we detail how the pointing model can be implemented. Finally, our 
work is summarized in Sec. 1 


2 Isotropic pointing model 


Basically, telescope mechanics are treated as a serial robot formed by two 
subsequent stages where the two stages are performing (nearly) orthogonal 
rotation. Throughout this paper we present the calculations in first equato¬ 
rial coordinate system, however, a similar approach can be involved for alt-az 
mechanics by simply replacing the variables appropriately. In the case of an 
equatorial mount, the first shaft (whose bearings are fixed to the ground) is 
responsible for the rotation around the hour axis while the second, perpendic¬ 
ular stage is responsible for setting the declination. In the following , we reca ll 
briefly the main ideas behind the isotropic approach used in Meszaros et al.1 
(1201411 . 


Throughout the derivation of our models, let us exploit the Cartesian rep¬ 
resentation of celestial coordinates. Let us define the standard direction of 
the telescope tube to the x+ axis, i.e. without any rotation, the telescope di¬ 
rected to the vector of po = (1,0,0). This is in accordance with the historic 
parameterization of the first equatorial system if we treat the hour angle r 
and declination S values as polar angles. However, one should keep in mind 
that this parameterization yields a left-handed (a;, y, z ) system of reference 
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(see Fig- HI) • An arbitrary point at the sky is then defined by the vector 


( cos S cos t\ 
cos <5 sin t 
sin 6 J 


( 1 ) 


while p(r = 0,5 = 0) = po- If one treats the whole telescope system as a series 
of active rotations, then one can write 


p(r, S) = P t • P d ■ po (2) 

where the rotation matrices Pt and Pd are defined as 

( cos r — sin r 0\ 
sinr cost 0 , 

0 0 l) 

( cos 5 0 — sin S\ 

0 10. 
sin S 0 cos 6 ) 

As it can be recognized, the assumption in Eq. @ is valid only in ideal circum¬ 
stances. In real applications, the various deflections presented in the telescope 
system are quantified by the rotation matrices H, X and T using the relation 

P '(t ,6) = (H • P t • X • P d • T) • po (5) 

Here, H quantifies the misalignments of the telescope hour axis with respect 
to the ground, X quantifies the misalignments of the declination axis with 
respect to the hour axis while T quantifies the misalignments of the detector 
or the tube with respect to the declination axis. These matrices of H, X 
and T are quite close to unity. Hence, a first-order series expansion can be 
involved in order to approximate these. The idea is based on the well-known 
relation between skew-symmetric and orthogonal matrices. Namely, for any 
skew-symmetric matrix A = —A 1 , the exponential of exp(A) is always an 
orthogonal matrix with unity determinant. In practice, let us write these in 
the form of 

/ 0 — c b \ / 1 — c b \ 

H = exp c 0 —a « c 1 —a I 

\—b a 0 J \—b a 1 J 

fO-fe\ ( 1 -/ e \ 

X = exp / 0 -d « / 1 -d] 

\-e d 0 ) y— e d 1 J 

/ 0 —ih\ / 1 — i h\ 

T = exp ( i 0 -g « \ i 1 -g\ 

\-h g 0 / \-h g 1 / 

We call the pointing model isotropic if the matrices H, X and T do not depend 
on the telescope position. 


( 6 ) 

( 7 ) 

( 8 ) 


( 3 ) 

( 4 ) 



4 



Fig. 1 The left-handed coordinate system implied by the first equatorial coordinates, i.e. 
the hour angle r and the declination 


Here we give the full expansion of P up to the first order. In the formula 
presented below, ct, St, c<j and Sd denote cost, suit, cos<5 and sin S, respec¬ 
tively: 


P 


( Ct^d St CtSd 

StQi StSd 

Sd o c d 


/ 0 0 0 

+ a I —Sd 0 Cd 

\StCd Ct St^d 


( Sd 0 Cd \ /-S t Cd —Ct s t s d \ 
0 0 0 I + c' I c t c d -s t —c t s d J 

-ctc d s t CtSdJ \ 0 0 0 / 


( s t s d 0 s t c d \ /c t s d 0 c t c d \ 

c t s d 0 — c t c d + e' s t s d 0 s t c d 
0 10/ \- c d 0 Sd / 


( 0 -CtSd St \ /-St -CtCd 0\ 
0 -s t s d -c t + i c t -s t c d 0 
0 c d 0 / \ 0 -s d 0/ 


(9) 


During the expansion, the matrices proportional to c and / are going to be 
exactly the same like the matrices proportional to e and h. Hence, the above 
equation depends only on c + / and e + h. Therefore, in the following we use 
c' = c + f and e' = e + h in our computations. 

If the expansion above (Eq. [9]) is multiplied by the reference vector po, the 
term proportional to g will also vanish. This term effects computations only 
when the information provided by the apparent field rotation is also included 
in the pointing model since the first column of the respective matrix is identical 
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to zero. All in all, the expansion of Eq. m is 

/c t c d \ / o \ / S d \ 

p'(t, S) = I s t c d + a I -Sd + b I 0 

\ s d / \C d St/ \-CdCtJ 



This equation can be written in the form of 

6 

P'O, S) = p(r, S) + OfcP k(r, 5), 

fc=l 


( 10 ) 


( 11 ) 


where the vector p(r, S) is defined by Eq. ([TJ) , the parameter vector of 
a = (ai,... ,ae) is equivalent to (a,b,c',d,e',i) and the pfc(r, 8) are defined 
accordingly to Eq. m- 

The parameters (a*,) can be obtained by a linear least squares regression 
if a series of observations is known with reported mount axis encoder posi¬ 
tions of ( T n ,S n ) as well as the respective apparent positions of (T n , D n ). Here, 
“apparent” refers to the hour angle and declination values after correcting for 
precession, nutation, aberration and refraction as well as one should accurately 
obtain the DUT1 = UT1 — UTC differen ces througho ut the derivation of the 
local sidereal time instances (see Wallace . 2002 , 20081 for more details). Once 
this series of 1 < n < N measurements is known, one has to minimize the 
merit function 

N 


X 


= £ 


Xn 


( 12 ) 


71—1 


where 


y 2 = 

A. n 


6 

£ 

k =1 


p(^n? d n ) T ^ ^ etfcPfc(T n , <5 n ) p(T n , D n ) 


Minimization of the above equation for % 2 yields a linear array of equations 
whose solution is then straightforward. Due to the three dimensional nature of 
the p vectors, this equation for \ 2 implies 3N independent terms. However, the 
effective number for degrees of freedom is only 2 N — 6 since for all 1 < k < 6 
values and for all observations, 


p{r n ,S n ) ■ Pfc(r„,<5„) = 0. (13) 

In other words, these 6 base vectors of p*, defined in the expansion m are 
always perpendicular to the observation direction p(r n ,<5 n ). 

In the following we extend the isotropic model described above with terms 
that provide significantly better accuracy. 
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Table 1 Coefficients of the expansion of the pointing model up to the first order in the 
spherical harmonics. Note that the terms marked with stars (★) and dots (•) are not linearly 
independent. Hence, the first order expansion of the pointing model has 22 independent 
coefficients instead of 24. 



3 Anisotropic extensions 


The idea behind the isotropic pointing model is the fact that the parameters 
ak of the respective expansion do not depend on the the pointing direction 
defined by r and/or <5. In order to attain better accuracy, we investigate now 
how the matrices H, X and T and hence the parameters cik can d epe nd on the 
position of the telescope. Recalling earlier works (see e.g. Zhang fc Wul . 12001 1. 
we examine the expansion of these terms via spherical harmonics. Namely, each 
of the Ofc constants are replaced by a linear combination of spherical harmonics 
up to a certain order L. Here L is the maximum order of the corresponding 
spherical harmonics and the respective Legendre polynomials (see also later 
on for the actual definitions). In practice, let us extend Eq. (fill) as 


6 l i 

p'(r,<5) = p(r,<5) + EE E ^k£ra^£rn (t ? 5)pk (t. S '), 

k= 1 £—0 m=—£ 


( 14 ) 
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where W m (r, 5) represent the real orthogonal spherical harmonics with the 
indices (£,to): 


Y / jx _ (Ke m P e m (sinc)) sin(mr) 
Cm ’ ^ Ki m P™ (sin 8) cos(tot) 

if to < 0 

if 0 > to. 

(15) 

In the above definition, the terms P™(-) denote the associated Legendre poly¬ 
nomials. The constants K( m can have arbitrary but non-zero values since in 
Eq. (1141). the Y( rn functions are multiplied by an unknown parameter. The 
series for the first three orders are 

loo = 1, 


(16) 

Yi,_i = cos 6 sin r, 


(17) 

Y\o = sin 5, 


(18) 

Yip = cos 6 cost, 


(19) 

Y 2j _2 = cos 2 8 sin(2r), 


(20) 

Y 2j _ 1 = cos 8 sin 8 sin r, 


(21) 

Y 20 = 3 sin 2 ,5- 1, 


(22) 

Y 2 i = cos 6 sin 8 cos r, 


(23) 

Y 22 = cos 2 5 cos(2t). 


(24) 


In Table [1] we also summarize the terms Ye m (T, 8)pk(r, 8) up to the order of 
t<L = 1. 


3.1 Parameter ambiguities 

It is essential to note that an important property of the terms Y( m {r, 5)pk (t, 8) 
is the ambiguity due to the lack of linear independence between some of the 
terms. For instance, up to the order of L = 1, there are two relations that 
should be taken into account: 


0 = P4 + P6*10, (25) 

0 = P 2 P 1.-1 + Pa^io + Pi^ii (26) 

These relations are valid for arbitrary values of (r, 8) and hence in the first 
case (Eq. l25l) either <2400 or agio should be omitted from the model fit. In the 
second case lEa. ECU) , one of the terms a 2 i t -i, 0310 or am should be omitted. 
We also mark these terms affected by the linear dependence in Table [TJ Due 
to the presence of these two equations, the L = 1 anisotropic pointing model 
has 6 ■ 4 — 2 = 22 free parameters in total. 



For completeness, we give here the six additional respective relations which 
appear if we consider spherical harmonics expansion up to the order of L = 2: 


P6f2,-1 + P4^1,—1, 

(27) 

P6f2l + P 4 Y 1 I, 

(28) 

P6^20 + P6 + 3p 4 Fi 0 , 

(29) 

3piF 2 l + p 3 + 3p2*2,~l + P.3>20, 

(30) 

3piF 22 - 6 pi + 3p 2 l2,-2 - 

— Pik 20 + 6 p 3 y 2 i, 

(31) 

3p 2 y 2 2 - 6p 2 + 3pil 2 -2 ^ 

-P2y>o + 6p 3 i2,-i 

(32) 


The omission of the respective aum parameters should be done according to 
these relations. For example, the first respective parameters in Eqs. (EZD m, 
i.e. a62,-i, 0621, 0620, ffli2i, o a i22 and 0222 can forcibly be set to zero during the 
fitting procedures (along with, for instance, the parameters O6io and am, see 
above). Therefore, the L = 2 anisotropic pointing model has 9-4 — 2 — 6 = 46 
free parameters in total. 


4 Interpretation of the model coefficients 

Our initial attempt was to create a pointing model that is completely based on 
a direct mathematical approach. However, the various quantities appearing in 
the expansion, i.e. the coefficients of akim as well as the respective functions 
of F^ m (r, 5)pfc(r, 6) can be assigned to various physical characteristics of the 
telescope mount. In the following, we briefly summarize these “roles” of the 
various akim coefficients. 


4.1 The isotropic case 


As it is known from earlier works (see e.g. Spillar et al. . 1993 : Buie . 2003lf 


the most essential parameters of a telescope pointing correction functions are 
related to the polar displacement, the encoder zero point offsets and the non- 
orthogonality deviation of the two axes. Indeed, it can be recognized that 
a = ai = ciioo and b = 02 = 0200 are the terms that are proportional to 
the polar displacement (in azimuthal and elevation directions, respectively), 
c' = <2,3 = 0300 is the hour angle encoder offset, d = 04 = 0400 is the deviation 
from the right angle between the two axes and e = 05 = 0500 is the declination 
encoder offset. The term i = ae = a 600 is the deviation of the optical axis in 
the direction of sidereal rotation. 
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4.2 Anisotropic models 

One of the most frequently cited anisotropic properties of a telescope mount 
arises from the gravity of Earth. Gravity yields a deflection of the telescope 
tube directly as well as indirectly via the unbalanced axes. As the barycenter 
of the various parts moves along as the telescope is slewing to various celestial 
positions, the net torque also results in anisotropic pointing offsets. In practice, 
the equatorial and declination stages are usually driven by worm gears. Hence, 
another kind of anisotropy appears in the mechanical system if we consider the 
ellipticity and/or eccentricity of the driven worm gears of these axes. Further 
details and examples about the interpretation of various anisotropic coefficients 
are found below, in Sec. 15.31 


4.3 Periodic errors 

An expansion by spherical harmonics can be used to include additional terms 
that have a dedicated physical relevance. A prominent example can be the 
characterization of periodic errors. Most of the classic telescope mounts are 
worm-geared systems and due to the limited resolution of the rotary encoders, 
raw hour angle and declination values are recovered from multi-turn encoders 
mounted to the worm shaft instead of the gear shaft. Let us suppose a gear 
ratio of G. As the term U 3 P 3 quantifies the encoder zero point offset in the 
hour angle, it can be deduced that the two additional terms of 

a 3,G,—GP3^G,—G + «3,G,GP3^G,G (33) 

characterize the periodic error. Namely, 

\J a 3,G,-G + a 3,G,G (34) 

is proportional to the amplitude of the periodic error while the angle 

arg(a 3 ,G,G,a 3 ,G,-G) (35) 

tells us the phase of the periodic error at r = 0 encoder position. 


5 Implementation and model tests 

In order to use the pointing model in practice, one should consider how it can 
be implemented for “real” telescopes. The actual implementation has to have 
two essential parts: 

— First, one should obtain an appropriately chosen series of astrometric cal¬ 
ibration frames of the given instrument in order to obtain the respective 
(T n ,S n ) and ( T ni D n ) values needed for the model. Afterwards, once these 
data are known, the parameters of the pointing model are needed to be 
fitted. 
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Fig. 2 The sky in the first equatorial system, as it is seen by the RCC telescope mount 
used for the evaluation of the pointing models. The hue levels (in parallel with the horizontal 
axis) represent the celestial hour angle values while the luminance levels correspond to the 
celestial declination values. The distance in the color space is equivalent to the apparent 
celestial distance while the projected distance in this graph is proportional to the movements 
of the equatorial and declination axes. The clear white stripe is the distorted projection of 
the north celestial pole while the fully black stripe (at the bottom and the top) is the 
distorted projection of the south celestial pole. The horizon is shown by the two parallel 
thick black curves while the various positive horizontal elevation levels (in steps of 10°) are 
shown as thin black dashed curves. The mount domain is limited by the thick red stripe. 
These limits correspond both to the horizon as well as the northern pillar of the equatorial 
axis pier. 


— Second, “implementation” also refers to the insertion of the model par am- 
eters into the telescope control system (TCS) in order to have a benefit 
during the observations. 

In this section, we summarize our method of implementation, considering both 
the model parameter regression as well as the integration in a TCS. 


5.1 Derivation of pointing model parameters 

To implement our model, one should employ a classic, purely linear least 
squares regression analysis on functions containing merely linear combinations 
of trigonometric expressions. By considering the definitions of the linear least 
squares method (y 2 , see Eq. [12]) as well as the actual model (see Eq. [TT]), one 
could derive that the merit function needed to be minimized is 

N 

x 2 = [p l=...(t„,<5„) - p(T n ,D n )] 2 , 

n— 1 


(36) 
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where 


P l=...(t,8) = p 


6 L 


J + XX X 

k —1 £—0 m=—£ 


^ktrnXim^J j <5)p fc (r,(5). (37) 


If L = 0, we get the isotropic pointing model (having 6 parameters in total) 
while for L = 1 and L = 2 , we get the anisotropic cases with 22 and 46 free 
parameters. Of course, the value of L can further be increased, however, one 
should keep in mind that additional identities similar to Eqs. (l?5l) - (E^l) and/or 
Eqs. (Em-® will appear during the expansions. 

Our current implementation for minimizing y 1 2 is based on shell scripts 
written in the language of bash. These scripts exploit standard UNIX text 
processing utilities as well as the If it utility of the FITSH package (I Pall . 
20121 ) . This implementation expects a four-column input file (where each line 
contains the T n ,S n ,T n , D n values) and it is available on the FITSH websit(0 
in the “Examples” sectioifl 


5.2 Integration in telescope control systems 

Once the pointing model parameters are known, they should be integrated in 
the TCS. However, such an integration needs a lower level access to the TCS 
software stack while the above described procedure for obtaining the coeffi¬ 
cients can be done by anyone who has observational access to the telescope. 
There are two ways to use the model: 

— First, Eq. (fTT!) (or Eg.[Tilwhen considering only anisotropic models) should 
be evaluated when one computes the apparent equatorial coordinates based 
on the encoder values. This type of calculation is performed when the TCS 
is asked to report the current position, i.e. after the reading of the encoders 
on the mount axes. 

— Second, the inverse form of Eq. (fill) needs to be evaluated during pointing 
and tracking of the telescopes. Namely, if one has to determine the en¬ 
coder values based on the apparent first equatorial coordinates. This type 
of calculation is performed when the TCS is commanded to slew the tele¬ 
scope to a certain position, i.e. before targeting the motors (and/or any 
mechanisms). 

Any software library that handles pointing model evaluation needs to be ready 
for implementing both of the above types. 

In practice, the inversion of Eq. (ED can be implemented by changing the 
sign in the series expansion coefficients. This type of inversion yields systematic 
offsets in the order of a kem . If the typical mount deflections are in the range of 
a few arcminutes (i.e. a milliradian), then the error caused by this improper 
inversion yields a systematic effect in the range of sub-arcsecond scale (i.e. 

1 http://fitsh.szofi.net/ 

2 http://fitsh.szofi.net/example/pointing/ 
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Fig. 3 Residuals of the pointing model applied to the Konkoly 1-meter RCC telescope. 
The left panel shows the respective components of the PL=o( r n, S n ) — p (T n ,D n ) vectors 
for the isotropic case while the right panel shows these difference vector components for the 
L = 1 anisotropic model. The scales for the difference vectors are the same for both panels: 
one degree in the main abscissa and ordinate refers to one arcsecond in the residual vector 
components. Note that the figures are redundant, not only because they show all of the 
(x,y), (y, z) and ( z,x ) component combinations but due to the fact that residual vectors 
are always tangential to the unit sphere. The strong structure of systematics in the residual 
in the isotropic case is quite explicit. Note also that the residual vectors would be as large 
as the individual panels if we plot the residuals without any pointing model. See text for 
further details. 


microradians) which is usually sufficient. If higher accuracy is needed, the 
inversion can be done in a two-step iteration, yielding errors in the range of 
a k£m’ Since practical applications do not rely on subsequent evaluation of the 
pointing model corrections, these errors do not accumulate. 


The above cited “Examples” section of the FITSH webpage contains a full 
implementation of an ANSI C library that can also be integrated in a TCS. 
In addition, this section of the webpage links some code snippets taken from 
the currently running version of an 1-m class telescope control software (see 
below the next subsection for more details). 
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5.3 Model tests 


In order to test the feasibility of our model we conducted a series of obser¬ 
vations with the recently refurbished 1-meter Ritcliey-Chretien-coude (RCC) 
telescope of the Konkoly Observatory. This telescope is installed on an English 
cross-axis mount. The northern pillar of this type of mount strictly limits the 
slew domain of the telescope within a certain range since the tube is unable 
to observe objects close to lower culmination. In the available telescope mo¬ 
tion domain we defined a grid of (T n ,S n ) values and exposures were taken on 
each grid point with a net time of 20 seconds. In total, we gathered 81 im¬ 
ages covering the sky nearly uniformly with the inclusion of polar crossing. 
Due to the partially cloudy weather, only 69 of these images have successfully 
been analyzed in order to have accurate astrometric solutions. The positions 
of these measurements are displayed in Fig. [2] Images were acquired with an 
Andor iXon-888 frame-transfer electron-multiplying CCD camera (EMCCD) 
and reduced_with the standard calibration procedures using the FITSH util¬ 
ities (fPdl l2012h . Since the net field-of-view of the optical setup is relatively 
small (3.4' x 3.4'), astrometric solutions were derived using images acquired 
in parallel with a guider telescope. For an initial sol ution, we ex ploited the 


offline version of the Astrometry.net package ( Lang et al. . 20ldh while the 


cross-matching of the two camera images as well as the final astrometric so¬ 
lutions were obtained with the approp riate tasks of the FITSH package using 


the USNO catalogue (IMonet et all 12003 1 as a reference. 


After obtaining the J2000.0 centroids of each frame, apparent first equato¬ 
rial coordinates (hour angle and declination) were obtained after correcting for 
precessi on, nu tati on, a berrat ion and refraction using the standard procedures 
( Meeusl . Il998t IWallaceL 12008 1. The list of encoder positions read at the expo¬ 
sure midpoint was used as an input series of (r n , S r ) values while the corrected 
astrometric solutions were used as an input series of ( T n ,D n ) values. 

The unbiased residuals of the isotropic, the first-order anisotropic and the 
second-order anisotropic pointing models were op = 69.1-10 -6 , <j\ = 14.8-10~ 6 
and (7 2 = 10.2 • 10~ 6 radians, respectively (where ox denotes the residuals 
corresponding to the Lth order in the spherical harmonics expansion). Since 
one arcsecond is equal to 4.84814 • 10 -6 radians, these values correspond to 
14.2", 3.0" and 2.1". We note here that without any pointing model applied 
(which is similar to setting forcibly all of the coefficients a,Mm to zero), the 
residual is op = 0.0011 radian, equivalent to 3.7'. This is ~ 16 times larger 
even than the residual corresponding to the isotropic (L = 0) case. We plot 
the typical structures of pointing model residuals for the L = 0 and L = 1 
case in the two panels of Fig. [3] 

The fitted values and the respective uncertainties of the model coefficients 
can also be used to measure the actual deflections of the telescope. For in¬ 
stance, in our case the isotropic model yields values for the polar displacement 
and encoder zero point offsets that can be used for tuning the telescope. The 
actual values for the polar displacement are a = 0.000163 ± 0.000012, b = 
0.000356±0.000010 while the encoder zero points are c' = 0.000844±0.000013 
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and e' = 0.000705 ± 0.000009 radians. It also turns out that the devia¬ 
tion of the two axes is also significant, being in the range of arcminute 
( d = 0.000315 ±0.000085). However, this value represents an average value for 
these deflections since the effects described in Sec. 14.21 can yield an anisotropy 
in any of the previously mentioned values. Indeed, for instance, the devia¬ 
tion of the two axes shows a clear dependence on the telescope position it¬ 
self, namely the corresponding coefficients are 0400 = 0.000573 ± 0.000091, 
a 410 = -0.000008 ± 0.000023, 041,-1 = -0.000077 ± 0.000056 and a 4 i,+i = 
0.000191 ± 0.000056 (also in radians). It means that the deflection of the two 
axes can vary with a full amplitude of at least 1.5 arcminutes on the total 
telescope motion domain. 

We also included two terms corresponding to the G = 240 gear ratio of the 
hour axis worm-drive of the RCC telescope. This inclusion yielded a decrement 
of the unbiased y 2 values of 3.5 (which should be compared to 2TV — P = 
116, where N = 69, the number of the frames and P = 22, the number of 
parameters in the L = 1 anisotropic model). At the first glance, this fact might 
prove that the pointing model is improved by this additional term. However, 
compared to the total degrees of freedom, this decrement is marginal as well as 
the respective amplitude that differs from zero only by 1.5-sigma. In addition, 
total measured periodic error of this telescope is < 0.15", which is well below 
the residuals of the model. All in all, we can conclude that this telescope cannot 
benefit with this inclusion, but other telescopes (e.g. amateur-class mechanics 
with typical periodic errors in the range of tens of arcseconds) might do so. 


5.4 Polar crossing 

It should be noted that spherical harmonics yield the same values while we 
observe the same celestial position before and after polar crossing. In other 
words, y^ m (r, 5) = Yi m (r + 180°, 180° — S) for all values of r and <5. However, 
some anisotropic physical effects (see Sec. 14.21 above) might have different 
yields depending on the pier side of the telescope. Since the telescope can reach 
most of the target positions in both configurations, it is worth separating the 
coefficients corresponding to these two configuration domains. Here we define 
two target points to be in the same configuration if the telescope can be re¬ 
positioned between the two points without polar crossing. 

The easiest way to do this is to consider points with 6 < 90° and 90° < <5 
in independent fits. Indeed, our analysis yields an unbiased residual of 1.6" 
and 1.7" for the L = 1 model if the two parts of the sky are treated separately. 
Both residuals are significantly smaller than the residual of 3.0" and even a 
bit smaller than the L = 2 model for the full coverage. If we perform this sep¬ 
aration in the isotropic model, the unbiased residuals are going to be roughly 
the same (namely, 14.0" and 12.1", in our case). 
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6 Summary 


In this paper we have presented a generic analytical anisotropic telescope 
pointing model derived from purely mathematical considerations. Our results 
show that this type of approach provides RMS residuals comparable and even 
better than earlier types of pointing models. It should be noted, however, that 
this pointing model does not account directly the temperat ure dependence of 


the telescope system. Based on earlier works ( Mittag et al. . 2008 1. the coeffi¬ 


cients for the drifts of the various parameters due to the thermal expansion can 
be as large as 0.2 — 1" per degrees Celsius. Hence, the accuracy of the ambient 
temperature should be within a few degrees Celsius in order to recover the 
pointing with an accuracy comparable to the RMS residuals. Uncertainties in 
the refraction model, for instance a variation of 15 hPa in the pressure and/or 
5 degrees Celsius in the temperature yield an offset of 1" — 2" even at moder¬ 
ate (25 — 45 degrees of) elevations above the horizon. Although an improper 
refraction model can be corrected via the Yg m functions, the correction fac¬ 
tors depend on the temperature and/or the density profile of the layers in the 
atmosphere. A major advantage of our model is the lack of parameterization 
via polar coordinates. In our tests, one of the grid points had a declination of 
</j 2000 = 89°55' and it is also a well behaved point, not an outlier - despite 
the fact that the respective value of 1/cos 5 is more than 700. In order to have 
a successful implementation for our model, one should employ only classic 
linear least squares regression analysis on functions containing purely trigono¬ 
metric expressions. The implementation of both the model fit and the model 
evaluation functions is straightforward, and is demonstrated on the FITSH 
webpage. 
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